Stability of Texture and Shape of Circular Domains of Langmuir Monolayers 
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Finite domains of a Langmuir monolayer in a phase with tilted molecules can be modeled by a 
simple elastic free energy of an XY order parameter with isotropic and anisotropic line tension 
terms. The domains can and often do contain nontrivial textures, which in turn influence the shape 
of the domains. Herein we investigate the properties of a simplified isotropic model with a single 
elastic constant. For circular domains a first-order phase transition is found between two distinct 
textures: an exterior defect (or "virtual boojum") texture, and an interior defect texture. Starting 
with a circular domain and either of these two textures as a ground state, we find shape instabilities 
develop that depend on the elastic constants and line tensions in the simplified model. In both cases 
a necessary but not sufficient condition for the onset of shape instabilities is the possibility for a 
local negative effective line tension to develop from the anisotropic line tension term. 

PACS numbers: 68.18.+P, 68.55.-a, 61.30.Cz, 61.30.Jf, 68.55.Ln, 68.60.-p 
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I. INTRODUCTION 

The biological importance of thin films of surfactants, 
such as DPPC (a primary component of human lung sur- 
factant p]]), and other potential commercial applications 
of self-assembled structures j| has spurred a resurgence 
of research in the study of Langmuir films. For the physi- 
cist, monolayers of surfactants at the air- water interface 
provide an interesting 2D system with very rich phase be- 
havior Techniques such as fluorescence microscopy 
and Brewster angle microscopy allow direct observation 
of the system as it proceeds through a phase transition. 
X-ray scattering can also be used. Atomic force mi- 
croscopy on Langmuir-Blodgett films can provide more 
detailed information that appears to be consistent with 
passive observations of the Langmuir film before extrac- 
tion §. 

Typically, as a film is compressed, it passes from a gas 
phase (G) to a more condensed liquid expanded phase 
(LE) to an even more ordered liquid condensed phase 
(LC) and then, sometimes, into a solid phase before fi- 
nally becoming a multi-layer film. Experiments have re- 
vealed a very complex variety of patterns, shapes, and 
textures present in the LC/LE and LE/G coexistence re- 
gions Upon rapid compression of the film, highly 
branched structures are often seen. These appear to be 
the result of diffusion limited aggregation H, and typi- 
cally they relax to more regular shapes. More modest 
compression rates tend to produce a finite number of 
condensed domains that grow in size, but not number 
||, as the film is compressed through a coexistence re- 
gion. Films composed of enantiomeric surfactants some- 
times exhibit domains with a preferred handedness (chi- 
rality) to their shape. The handedness alternates with 
the choice of enantiomcr and in racemic mixtures the do- 
mains are achiral (|Jill ■ Within domains of the ordered 
phase, complex structures can sometimes be observed as 
well • The addition of cholesterol to these systems 
often leads to labyrinthine or striped patterns, and to 



spiral domains of a preferred handedness Pj . 

Molecules in LC phases of Langmuir films either align 
normal to the air-water interface, or they can tilt rel- 
ative to the interface ]F2||l5| ], thereby defining a two- 
dimensional vector in the plane of the film. In addi- 
tion some films exhibit hexatic order as well as tilt order 
]l4] , [i"5| ]. Domains of both tilted and non-tilted phases 
exhibit non-circular shapes P,pdpp|,p^| . Static and dy- 
namic properties of non-tilted domains as well as mod- 
ulated equilibrium phases are well explained by dipole 
interactions among the electric dipoles aligned normal to 
the layers |]l7|-^]. Domain shapes and textures in tilted 
domains clearly depend on the existence of orientational 
order in the plane of the film (l^j2^]. For example, it 
is difficult to explain chiral shapes that appear in chiral 
films by a dipole model. 

In this paper, we will consider a particularly simple 
model for tilted domains in which tilt order is described 
by an isotropic XY-raodel and in which coupling between 
texture and domain shape arises because of an interaction 
favoring alignment of XY- vectors at a specific angle with 
respect to the domain boundaries. This model is clearly 
oversimplified, but it can explain observed shapes and 
textures in some systems . A complete model 

would allow for different elastic constants for splay and 
bend distortions and for splay-bend coupling in chiral 
systems pq| in addition to interactions between electric 
dipoles with components both parallel and perpendicular 
to the film. Nonetheless, this simplified isotropic model 
predicts non-trivial shapes and textures, the full range of 
which have yet to be analyzed. In spite of its simplicity, it 
involves highly non-trivial couplings between texture and 
shape that need to be understood before more realistic 
models can be studied. 

Our model is characterized by three energy parame- 
ters; the isotropic elastic constant K (units of energy), 
the isotropic line tension 7 (units of energy/length) and 
the line tension /i favoring alignment of the XY- vector m 
with the tangent to the domain perimeter t. We obtain 
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the following results for this model. For /z/7 < 1, circular 
domains are globally stable. Circular domains of radius 
R can have two equilibrium textures: one with a strength 
+1 disclination [E7J with core radius £ at the center of 
the domain and one described by a +2 defect (virtual 
boojum) exterior to the domain [^8[ . The interior de- 
fect texture is favored for large values of R/£ whereas 
the exterior defect texture is favored for large values of 
K//j,^. For a circular domain, the exterior defect texture 
never becomes locally unstable. The interior defect tex- 
ture does, however, become locally unstable with respect 
to moving the defect toward the edge of the circle. When 
/1/7 > 1, the circular shape and associated textures can 
become unstable. When the defect is in the interior, the 
circle can become unstable with respect to n-fold mod- 
ulations as the domain radius R is increased. The value 
of n that first becomes unstable depends on /i/7. When 
the defect is outside, the domain becomes unstable with 
respect to n — 2 distortions first, as R is increased. 

We will proceed with our investigations as follows. We 
will review the interior and exterior defect textures for a 
circular domain and then, comparing them, will present 
a phase diagram for these two textures. Next, we will 
examine the shape stability by taking an initial config- 
uration consisting of a circular domain and one of the 
two extremal textures. For each of the two textures we 
will then allow for perturbations in the shape and tex- 
ture, calculate an effective free energy dependent only 
upon the shape, and examine the stability of the circular 
domain with respect to shape deformations. 



II. THE MODEL AND OUR AIMS 

Henceforth we will be concerned only with the follow- 
ing problem. Find the minimum of the free energy p9[ 
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K, 



K b (V x rh') 2 



d 2 x 



—fj/ d> rh' • ids — if &> rh' x ids + 7 <z> ds, (1) 

JdD JdD JdD 

over all domains D with fixed area A (= irR 2 ), and over 
all unit vector fields rh', for given values of the elastic 
constants (K s , Kb, fi', 77' , 7). K s and Kb are the usual 
bend and splay elastic constants, 7 is a line tension, t is 
the tangent vector to the curve dD, and jjl and rf are the 
coefficients of spontaneous bend and splay, respectively. 
We have chosen to write the spontaneous bend and splay 
contributions as line integrals using 



V • m'd 2 x = d> m! x ids 

D JdD 

V x m'd 2 x = <p rh' • ids. 

I • JdD 



(2) 



This allows us to identify — //rh' ■ t — 77'm' x t + 7 as an 
effective anisotropic line tension. If 7 < y /1' 2 + rf 2 , then 



choosing rh' to follow t for a given curve dD will allow 
the anisotropic line tension to be negative. We will see 
that within the context of the one-coupling-constant ap- 
proximation [K = K s = Kb), allowing the effective line 
tension to become negative will be a necessary condition 
for instabilities in the shape of a circular domain. 

We will now restrict our attention to the one-coupling- 
constant approximation, simplifying the form of the bulk 
free energy. Defining $' via 



rh' = cos $'e x + sin &'e y , 



(3) 



we have 



D 

K 



K s (V • rh') 2 + K b (V x rh') 2 



V$' ■ y^>'d 2 x. 



d 2 x 



(4) 



In appendix |b] we show that the transformation 



$= $' - tan 



-1 V_ 
V' 



TV 
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M= V V 2 + V 2 

rh= cos $e T + sin <f>e„ 



(5) 



simplifies the free energy further: 
F=-K [ {\7&) 2 d 2 x -[il rh • Ids- 

2 J D JdD 



+ 7 <z> ds. (6) 

JdD 



We could just as easily have transformed away the spon- 
taneous bend term in favor of an effective coefficient for 
spontaneous splay. Thus, although Eq.(|6|) is a chiral free 
energy (the spontaneous bend term is chiral), chirality 
can only manifest itself through the texture rh and not 
through the domain shape D {D is invariant under the 
transformation, see App. |b] ). 

To make contact with previous work |T(| we note that 
we could have written our surface term in the form 



/ cr(0 - <$>)ds, 
JdD 



(7) 



where the outward normal to the boundary is given by 



n = cosfe T + sin fa 



and 



<t{x) = 7 + (a„ cos nx + b n sin 1 



(8) 



(9) 



Taking b\ — —\i with all other a n 's and b„s set to zero 
will recover our effective line tension 



7 ® ds 

IdD 



H <x> rh • ids. 

IdD 



(10) 



Similarly a nonzero a\ will produce a spontaneous splay 
contribution. The transformation (0), in this language, 
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allows the elimination of the a\ (b±) term in favor of a 
new effective b[ (a^). 

Rudnick and Bruinsma |l6| have demonstrated that 
a nonzero ct2, in addition to a nonzero a±, leads to 
non-circular shapes and nontrivial textures in agreement 
with some experiments pBf . However, Riviere and Me- 
unier [ [Tl| can account for these same shapes and tex- 
tures without a nonzero 02 by incorporating the appar- 
ent anisotropy between K s and Kb- We will show that 
even with K s = Kb and 02 = [as in Eq. (^)] shape in- 
stabilities still exist when 7 < fj,, that is, once a negative 
effective line tension is allowed. 

To see that instabilities should arise, we first identify 



A,, 



K 
~j~iR 



7 



(11) 



as two independent dimensionless parameters. Ab is a 
measure of the competition between the bulk energy and 
the anisotropic piece of the line tension. For a given dD, 
Af, will determine the preferred texture. A s measures the 
relative strength of the anisotropic and isotropic contri- 
butions to the effective line tension. When A s > 0, the 
effective line tension can be negative. 

When Ab — > 00 we can assume that K — > 00 and hence 
that m will be constrained to be a constant. Thus, with- 
out loss of generality, we can assume that rh = e,,. This 
form for rh yields 



rh • ids = 



OD 



(12) 



for all dD, and as such the value of \x is irrelevant. In 
fact, the free energy is simply 



F = 7 i> ds. 

'dD 



(13) 



For fixed area, A, the preferred shape will simply be a 
circle. 

In contrast, when Ab — * we assume K — > and thus 
rh can take any form it chooses in the bulk. In particular, 
we can always choose rh such that rh • t = 1 at all points 
on the curve, for any curve dD. Thus now our free energy 
becomes 



F = (7 — fj) <b ds. 

'dD 



(14) 



If 7 > n then certainly we again have a circle as the pre- 
ferred shape. But if 7 < /i then the domain attempts to 
maximize the length of its boundary, dD. This results in 
pathologically distorted shapes, or unbounded domains. 
Thus, when 7 < /a we anticipate that for large Ab the do- 
main shapes will be circular, but for small Ab the shapes 
will begin to distort. However, when 7 > ji it would not 
be at all surprising to hnd that the circular domains are 
preferred for all values of Ab. 



We are now ready to begin our search for optimal shape 
and texture configurations. Although we will restrict 
our attention to the one-coupling-constant approxima- 
tion, which provides us with a particularly simple bulk 
Euler-Lagrange equation 



V 2 $ = 0, 



(15) 



it is still not a trivial task to find the optimum $ for a 
generic domain D. We will approach the problem as fol- 
lows. Noting that experimentally the domains appear to 
be circular for small A, we will first fix the domain shape 
to be a circle and try to find the optimum 4>. This is 
tantamount to assuming that 7 is infinite. Having found 
the optimal texture for a circular domain, we will then 
allow the shape and texture to deform. The next section 
will deal with the task of finding the optimal texture for 
a circular domain and the succeeding section will deal 
with instabilities in the shape of the domain. 



III. THE OPTIMAL TEXTURE FOR A 
CIRCULAR DOMAIN 

The main question to address here is whether the do- 
main should contain a topological defect. In experiments 
by Riviere and Meunier [[llj the domain clearly does not 
contain a defect and has the "virtual boojum" texture 
previously described by Sethna |^8| . However in (!(]] Qiu 
et al. observed a star defect inside the domain, which 
indicates that it may be possible to have an isolated +1 
defect inside a domain when hcxatic order is not present 
(as the hexatic order is lost the arms of the star defect 
retract and one is left with a simple +1 defect pp[-p2|). 
Also, in thin films of smectic-C liquid crystals, which can 
be approximately modeled by the same free energy, iso- 
lated + 1 defects are often observed Jj~7| . 

Fixing the domain to be a disk, the isotropic line ten- 
sion term will be irrelevant to this discussion since it is 
unaffected by changes in texture alone. The anisotropic 
piece of the boundary energy, however, favors having rh 
parallel to the tangent to the boundary and thus prefers 
having either a +1 defect inside the domain or a +2 defect 
on the boundary (sometimes referred to as a boojum). 
The bulk portion of the free energy clearly prefers that 
rh simply be uniform within the domain. Thus, whether 
it is preferable to have a defect inside the domain or not 
depends upon the competition between the bulk elastic 
energy and the anisotropic line tension, as well as upon 
the energy cost required to nucleate a melted defect core. 

In this section we will first review the virtual boojum 
texture, which is an extremal texture for the case where 
ther e is n o defect inside the disk (l6). In fact, in sec- 
tion IV A we will show that this is a locally minimal 
texture with respect to our free energy. Next we will 
consider textures associated with the defect inside the 
domain. We will note that the texture for an isolated +1 
defect at the center of the domain is extremal and for 
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A\, < 1 it is a local minimum as well. We will investi- 
gate a class of ansatz textures consisting of an isolated 
+1 defect located anywhere within the domain. We will 
find that within this restricted class of textures, there is a 
first-order transition between the texture with the defect 
at the center of the disk and a texture with the defect 
close to the boundary. Thus, we find that not only is the 
texture of the defect at the center of the disk unstable 
for Af, > 1 but that it also fails to be a global minimum 
for even smaller values of A b within the class of textures 
with an enclosed defect. In contrast, the virtual boojum 
texture is a local minimum for all values of A b . 

Finally, we compare the virtual boojum texture with 
the texture of the defect at the center of the disk. We find 
a first-order transition between these two textures which 
is most easily described by the dimensionless parameters 
c = K/2[i£, x = R/£, and g, where £ is the coherence 
length and g is a phenomenological parameter associated 
with the energy cost of a defect core (in our calculations 
we will take g = 1). 



A. Virtual Boojum 
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(a) (b) 
FIG. 1. Textures arising from a single +2 defect. The con- 
stant term in <E> has been chosen to make rh as close to t (the 
tangent vector) as possible. With the defect at the boundary 
(a) we have rh • t = 1 but we also have very large gradi- 
ents in rh near the defect. In (b) these large gradients have 
been expelled along with the defect but now we no longer 
have rh ■ t = 1. The competition between the ideal boundary 
(rht = 1) and bulk (dim = 0) conditions leads to a preferred 
separation distance between the domain and the defect. 

Let us now recall the virtual boojum texture ^8|. This 
is simply the texture produced by a single +2 defect sit- 
ting outside the disk (Fig. [I]). With the following nota- 
tion, 



the field for a +2 defect located 
Yd = {r d cos (j> d ,r d sin <j) d ) = (x d ,y d ) is just 



$ = 20 rd (r)+C. 



at 



(17) 



Clearly the texture satisfies the bulk Euler-Lagrange 
equation, V 2 $ = 0. Furthermore, this collection of tex- 
tures (one for each and C) smoothly varies in its ex- 
tremes from the uniform texture (r^ — > oo, rh — > con- 
stant), which is preferred when \x — 0, to a texture where 
rh • t = 1 (r d = R), which is preferred when /i — > oo. 

Taking the domain to be a disk of radius R centered 
at the origin we then find that Eq. (||) yields 



2-k—R 2 sin(C + 4> d ) 
rd 




+2irjR, 



(18) 



where we have assumed that the defect is outside the 
disk at least a core radius, £, away from the boundary 
(rd > R + £). Minimizing over C we find the preferred 
value of C to be, 



(19) 



Without loss of generality we can take 4> d = 0, and finally 
write the energy for a +2 defect outside the disk, on the 
positive x— axis as, 



2ir—R 2 + 2irjR. 

rd 




For a given value of the dimensionless quantity 

K 

A b = — 
fj,R 

there is a preferred stable minimum value 
given by 



r° d = R(A b + Jl + Al) 



(20) 

(21) 
,Hl for r d 

(22) 



Furthermore, we see that using this value of Td in our F + 2 
above yields F +2 (r°) < 2ir-fR, revealing that this texture 
is preferred over the uniform texture for all values of A b . 

Rudnick and Bruinsma |l(| have shown that this tex- 
ture with the defect at r° is in fact extremal. We will 
further show that this texture is indeed a local minimum 



when we investigate shape instabilities in section IV 



r = (r cos 4> 1 r sin <f)) — (x, y) 
z = x + iy 



(16) 



'r,, 



(r) = Im(log \z — Zd\) = tan 



Vd 

Xd 



4 




r d = (0,0) r d = (0.75 R,0) r d = (0.75 R,0) 

(a) (b) (c) 

FIG. 2. Texture arising from a +1 defect inside the disk. The textures in (a) and (b) arise solely from the interior defect, 
whereas the texture in (c) includes the effects of an image defect outside the domain. Notice that in (a) and (b) we have 
m ■ t = 1 whereas in (b) we do not. 



B. Defect Inside 

Here we will look at the class of ansatz textures 

7T 



(23) 



where we assume rd < R — £ (i.e. the defect is inside the 
disk and at least a core radius, £, away from the bound- 
ary). Without loss of generality we take 4>d = and ac- 
cordingly we have already chosen the constant term to be 
the optimal value. Again, clearly, all such textures satisfy 
the bulk Euler-Lagrange equation, V 2< £ = 0, within the 
domain, except at the position of the defect. We note 
that when = 0, m is parallel to the tangent to the 
boundary (Fig. ^). 

Applying Eq. (||) we find that the energy of this tex- 
ture is 



F+i — ttK log 
+2njR - 



R 



log 1 



' d 
R 2 



4*1! (2) 

(24) 



where we have introduced a term for the core energy, 
e corc . E(x) is the complete elliptic integral of the second 
kind, 



E(x) 



(25) 



We will see in section [V , when we investigate changes in 
shape and texture, that when r c i — and A b < 1 this tex- 
ture is a local minimum with respect to all variations in 
the texture (not just variations within this ansatz class). 
However, we find that rd = and A b < 1 does not always 
yield a global minimum even within this class (p3L). 



The free energy ([24|) has an interesting behavior as we 
vary rd- As we move the defect off center, there is an 
energy increase from the boundary term but a decrease 
in energy from the bulk term. The two contributions 
have quite different functional forms, and we find that, 
depending upon the value of At, = K//iR, it may be 
preferable either to keep the defect at the origin or to 
allow it to migrate towards the boundary. Maintaining 
< rd < R — £ we find the following behavior of fl24| ) 
within this class of textures, 



A b < 0.23 
A b < 1 
A b > 1 



rd = is a global minimum 
rd = is a local minimum 
rd = is a local maximum. (26) 



This behavior is summarized in Fig. ||. We see that 
within this class of textures there is a first-order tran- 
sition in the position of the defect as we change A b . 
When rd = becomes the global minimum, it may in 
fact already be the case that the energy barrier that must 
be overcome in order for the defect to migrate from the 
boundary to the origin may already be too high for this 
to be a physically probable event. Note that we have lim- 
ited our attention only to a particular subset of potential 
textures and this barrier may not remain if we were to 
consider all possible textures. However, we are aware of 
no experimental observations corresponding to a single 
+1 defect near the boundary of the domain. Thus, over- 
all we conclude that when a defect is present within the 
domain it is very likely to be at the origin in a stable 
local minimum. 
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FIG. 3. Plots of F+i/K vs r d /R for a disk of radius R for 
several values of At- We have taken £/R = 1/1000 and have 
restricted rd/R to be less than 1 — 1/1000, so that the en- 
tire defect core is within the disk. We have also not included 
the core energy e core or the isotropic line tension 2-k^[iR in 
F+i, which here are merely uninteresting constants. Note 
that Vd = is a local minimum whenever At < 1 and ~ R 
appears to be a local minimum for all of the plots. When 
At < 0.23 then rd — appears to be the true minimum. 
Notice that the barrier separating the two minima is approx- 
imately 4K at A b = 0.23. 



4/R=l/1000 




FIG. 4. A comparison of the energy of an isolated +1 defect 
in the disk, F+i with the energy of a +1 defect inside along 
with an image defect outside enforcing strong pinning bound- 
ary conditions (rht = 1), F+i.+i. Note that the energies are 
equal at R d = 0. At the boundary of the disk (rd = R — £) 
_F+i is always lower. Even though for some intermediate val- 
ues of r d f+1,+1 is smaller than F+i it is never substantially 
smaller. 



We note that one can also consider the texture, 



4> 



0r d (r) +c/) rdR 2 /r 2(r) 



2' 



(27) 



which corresponds to a defect inside the domain along 
with an image defect outside the domain which enforces 
strong pinning boundary conditions, m • t = 1 (where 
again we take 4>d = 0). Within this class textures the 
defect always prefers to be at the origin, and when Ab > 
0.23 this texture provides at best a modest improvement 



over (|23j) and never offers a lower global minimum (Fig. 

Of course, when = this texture is the same as 
(|23|) with rd — 0. 



C. Which texture wins? 

We have presented two textures which are local minima 
of the free energy (0). For the +2 exterior defect texture 
we found a local minimum for all values of A\ ( ( |r7j ) with 
4>d = and Td = r~d ( |22] ) ). For the +1 interior defect 
texture we found a local minimum only for Ab < 1 ((|23|) 
with rd = 0). Furthermore, we found that this local 
minimum was not the global minimum even within the 
class of ansatz textures (^3|). Here we will compare the 
energies of the two locally minimal textures to see which 
is favored. 

For small R the bulk energy dominates the boundary 
energy and we expect the uniform texture to win. Since 
the +2 defect texture produces the uniform texture in 
the limit Ab — ► 00 (i? — » 0) we anticipate that for small 
R this texture will be favored. For large R (Ab — > 0) 
both textures satisfy strong pinning boundary conditions 
and thus have equal contributions from the surface en- 
ergy. However, the bulk energy cost of the exterior defect 
texture is greater that that of the interior defect texture. 
This can be understood qualitatively by recalling that the 
bulk energy cost of a defect of charge q centered in a disk 
of radius R is nKq 2 \og(R/£) (§7|. In tne umit of l ar S e R 
the +1 defect energy will be TiK\og(R/ whereas the +2 
defect energy will be w (l/2)2 2 7riClog(i?/C). The factor 
of (1/2) arises because the +2 defect is at the boundary 
of the disk, not the center. 

Thus, we will find a first-order transition between the 
two textures. In discussing this transition we find it more 
convenient to refer to the dimensionless quantities 



K 



and 



R 

X = I 



(28) 



rather than A^ and R/£. Assuming that R > £ we find 
that when c < 0.85 the interior defect solution is favored. 
For larger values of c, there is a critical radius, i? C rit(c), 
above which the interior defect texture is favored and 
below which the exterior texture is favored. Thus i? C rit(c) 
defines the first-order phase boundary T in Fig. |[ 

Note in Fig. || that we have also included the curve 
Ab = 1 which is the boundary at which the interior de- 
fect solution becomes unstable. There is a large separa- 
tion between this boundary and the transition line T. In 
the region above the At> = 1 curve both textures are local 
minima and hence a physical system may well find itself 
trapped in a meta-stablc local minimum when Ab < 1. 
Furthermore, in section III B we found some circumstan- 



tial evidence suggesting that the barrier preventing the 
expulsion of the interior defect may be quite large for 
physically realisable values of Ab (Fig. |3|). Indeed in an 
experiment |25[| that apparently has the exterior defect 



G 



texture, taking £ « lOOnm yields c « 2. However, the 
radii of the domains appears to be as large as 10/im, 
which is much larger than the critical radius for c = 2 
(-Rcrit(2) « 2pm), suggesting that it may be the case that 
this texture is a meta-stable equilibrium state. 




FIG. 5. The two figures display at different scales the phase 
diagram for the +1 defect at the center of the disk versus 
the exterior +2 defect texture. The curve T is R = R C rit(c) 
where i? cr it (c) is the c ritical radius for a given value of c as 
discussed section [II C . The lower curve is At, = 1. To the left 



and above T the interior defect texture is favored, below the 
exterior defect is favored. However, the interior defect texture 
is a local minimum in the region above the curve At, = 1 and 
the exterior region is a local minimum throughout the region. 

To identify the transition we simply compare F+i and 
Taking e corc = nK and noting that A\ = 2c/ x we 
can write (EHI) as 



2c 2c 
F + i = irpR, I — log x - 2 H 

X X 



2-k^R. 



(29) 



With 



and Td 



we can also write (EOh as 



F 



+ 2 



where 



TTflR 



4c. 



log 1 



2h 



2c 



h(y) = 



y + V 1 + y 2 



(30) 



(31) 



Examining F + i — F + 2 we find that when c < c* rj 0.85 
the interior defect solution is preferred for all values of 
R (or equivalently all values of x since we are assuming 
that £ is constant). For larger values of c we find that 
there is a critical value of R, i? cr it(c), below which the 
exterior defect solution is favored (Fig. 



IV. SHAPE CHANGES 

For circular domains we have found two distinct tex- 
tures which are local minima of the free energy (this will 
be demonstrated soon). Now for each of these textures 
we would like to know if the configuration (D, $ )j where 



D is a disk and <!>o is either of the two textures, is a min- 
imum with respect to changes in shape and texture. We 
will find in each case that although the texture is sta- 
ble when the shape is fixed, the configuration (shape and 
texture) is unstable for certain values of At, and A s . We 
will see that when p, > 7 (A s > 0) instabilities will begin 
to arise. 

To examine the stability of the configurations we need 
to introduce perturbations in the shape and texture. To 
this end we take D to be the region bounded by the curve, 
given in polar coordinates by 

r(<t>)=Ro+ ]T Pn e in * 

neZ-{0} 

= Ro + (a ra cos ruj) + b n sin rwj)), (32) 



neZ+ 



where 



Pn = P*- n = ^(a n - ib n ) 



(33) 



and R is defined by 



R = R\1- Yl 



PnP- 



neZ-{0} 



R 2 



(34) 



al + hl 



2 ^ R 2 

neZ+ 



With this parameterization the area of the domain is sim- 
ply 



A = ttR 



(35) 



Thus keeping R fixed keeps the area constant. Note that 
we have included the n = 1 terms in (jH). Although 
these modes can usually be ignored as being just trans- 
lations of D (to lowest order in the Fourier coefficients), 
here the texture distinguishes a particular origin for the 
coordinate system and as such translations of D are as- 
sociated with nontrivial changes in the configuration (see 
Ap P . 0). 

To allow for perturbations in the texture we will take 



$(r) = $ (r) + Y r l " l P„e in0 with P n = PI 



neZ-{0} 



(36) 



where $o( r ) is one of the two textures found in section III 
( @ and @) and the perturbation is the most gen- 
eral nonsingular addition allowed by the Euler-Lagrange 
equation V 2 $ = 0. Taking 0{P n ) = 0{p n ) = 0(e) we 
will proceed to calculate F to order e 2 . This will provide 
us with an approximation to 



F(p,P,K, 7 ,p,R), 



(37) 
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where 



P = (pl,P-l,P2,P-2,---) and 
P = (P 1 ,P_ 1 ,P 2 ,P_ 2 ,...). 



(38) 



Minimizing this over P we will obtain the effective free 
energy, 



F cS (p, K,j,fi,R), 



(39) 



as a function of the elastic constants, the radius, and the 
shape perturbation parameters correct to 0(e 2 ). 

In both cases we find that F e g contains no terms of 
0(e), that is no terms linear in p n , demonstrating that 
the configurations are extremal. We also find that when 
fj, > 7 shape instabilities will be present for certain values 
of A\, and A s . In fact, for any given A s > we find that 
instabilities arise for sufficiently small A^s. 

Ideally we would like to calculate F e ff correct to all 
orders in {p n }- Then minimization of P e ff over the p n 's 
would yield distorted equilibrium shapes. Unfortunately 
this is a prohibitively difficult task. Alternatively, we 
would like to calculate F c g to high enough order in e such 
that the effective energy is bounded below. This would 
provide approximations to the equilibrium values of the 
Pn's, and thus approximations to distorted equilibrium 
domain shapes. Calculating F c g to 0(e 4 ), while difficult 
but not impossible, is not sufficient. This appears to be 
due to the form of 7 § dD ds as a function of the p n 's. 
Although this isotropic line tension term should be sta- 
bilizing the boundary against deformations, the fourth 
order coefficients are negative and as such they lead to 
F e ff being unbounded below at 0(e 4 ). Because of exten- 
sive mode coupling and the exponential increase in the 
number of terms in F e g we have not proceeded to 0(e 6 ). 



Noting the sin(C + 4>d) term in (y_8J), we see that ro- 
tating the defect about the center of the disk and simul- 
taneously adjusting C to maintain C + <j)d — —n/2 costs 
no energy. As we have just noted, an appropriate choice 
of the P ra 's will produce this rotation, and thus if C were 
allowed to vary, we would discover this Goldstone mode 
eventually. By fixing C to be — ir/2 we freeze out this 
mode. Experimental evidence suggests that this is the 
appropriate course to take. The texture associated with a 
+2 virtual defect has been observed ||n] , ^5fl , demonstrat- 
ing that on the time scale of the observations the virtual 
defect does not freely rotate about the domain. If the 
defect did execute such a motion, then the time averaged 
texture would appear as the trivial texture. Though a 
careful consideration of the dynamics might be in order, 
in this paper we will not concern ourselves further with 
the origin of the observed angular stability of the virtual 
defect. 

We find it more convenient now to use the real rather 
than the complex fourier coefficients in calculating F: 

A n = 2Re[P n ] B n = 2Im[P„] P„ = ~(A n + iB n ) 

a n = 2Re[p„] b n = 2Im[p n ] p n = -(a n + ib n ). (42) 

Thus, now A = (A 1} A 2 , . . .) and B = (B 1 ,B 2 , . . .) 
describe texture changes whereas a = (ai,a 2 , . . .) and 
b = (61, b 2 , . . .) describe shape changes. Taking K = 1 
and R = 1 to set our energy and length scales, and per- 
forming a plethora of contour integrals we eventually find 

P(A, B, a, b, A b , A s ) = {A\H AA \A) + (B\H BB \B) 
+ (a|P aa |a) + (b|P bb |b) 
+ {A\U Ab \b) + (B\U Ba \a) (43) 
+ F, 9 , 



A. The exterior defect 



where 



We will now apply the outlined procedure to the case 
of the exterior defect. We take D as in ( |32"|) and we take 



$ o = 20 r o(r)--, 



(40) 



the preferred texture for a disk, as found in section III A . 
We will not allow r° to vary, nor will we allow for varia- 
tions in the constant term (C = — 7r/2). Allowing to 
vary would present us with a redundancy in the definition 
of $ (|3^). Namely, changes in r° are already accounted 
for by the P„ 's in (B6h . If we wished to have 



(41) 



where ^ this could be achieved by an appropriate 
choice of the P„'s in ( |36| ) since the perturbative terms 
form a complete set of solutions to Laplace's equation in 
the disk. 



(HAA)n,m — (Pss)n,m — 

— pm ( I 1 



(44) 



4 V (rjj) m + n-1 ( r ^m+n+l 



4 VW) 
■ 2r 0°m,n- 



0\\m—n\—l 

2Trm6 m 



jyoym-nl+i 



(UAb)n,m — 
pTT I 1 



n—m+ln 1 ^' 



(UBa)n,r, 



pYK 



\ ^O^m+n 



(45) 
(46) 

, / pn-mn 1 ~ Sm > n + OmS 1 



and 
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(Haa)n,m = -(h(n,m) + h(n, —m)) + -jK-f{n 2 - l)6 m ,n, 

(47) 



(Ha,) 



i,.m = -j{-h{n,m) + h{n, 



-TO)) + -7T7(ri 2 



where 

h(n,m) = 7174(1 - <5 m _„) 
7i + n| — 2 



l)<5m,n, 

(48) 
(49) 



|to + n| — 2)(|m + n| — 1) \m + n\(\m + n\ + 1) 



2 ( r O)|m+n|-l 



2 ( r 0)|m+n|+l 



2tt 



{r° d ? 1 
2tt 



|m + n\ + 1 



( r 0)|™+»| ^ (r o )a _! ((rO)a-i)5 

For convenience we have chosen to leave the above in 
terms of /i, 7 and r°. With if = 1 and R — 1 we have, 



1 



A 2 



1 = 



-1 



(50) 



A s + 1 A s 



1 



Note the absence of any terms linear in {A, B, a, b} in 
( p3| ) . This confirms that the unperturbed configuration is 
an extremum. Furthermore, Haa, Hbb, H aa and H bb are 
all positive definite. Thus the unperturbed configuration 
is actually a minimum with respect to changes in either 
texture or shape alone. Any instabilities must therefore 
arise as a result of the couplings, U Ab and Us a , between 
shape and texture. 

Before minimizing over A and B to find F e g we note 
that two Goldstone modes are still present. A simulta- 
neous translation of the domain and the defect will not 
change the energy, F. This implies the existence of 2 zero 
eigenvalues in F, one for horizontal and one for vertical 
translations. These translations are achieved through ap- 
propriate choices for A, B,a and b. For example, the 
infinitesimal translation ee x is achieved by taking 



2e 



(»■§) 



0\n+l 1 



(51) 



The two Goldstone modes should be present in F e s and 
we will use their existence as a check on our calculations. 
Symbolically, calculating F e g is straightforward, 



where 



- U Ba H BB U Ba 



Qbb = ~~~:UXbH A ^UAb- 



(53) 



However, the actual calculation of F e s is not so easy. 
Neither Haa nor Hb b is diagonal, the off-diagonal terms 
are not in general small, and the higher order harmonics 
(in A and B) are not necessarily unimportant. We will 
resort to numerical procedures to aid us in finding the 
eigenvalues and eigenvectors of -F e ff- 

If r° d is large then the higher harmonics in A and B do 
become less important, and we anticipate that working 
with finite dimensional Haa's and i/ss's will provide a 
reasonable approximation to F c q. Furthermore, if 7 ^ 
then the higher harmonics in a and b are also rapidly 
suppressed in H aa and H^, allowing us to work with 
finite dimensional 7J aa 's, H^s and U's as well. 

Noting that a = (1, 0) and b = (1, 0) are the Goldstone 
modes for F c ff provides us with one gauge of how well we 
have approximated F c s with finite dimensional matrices. 
Even if we do include too few modes to correctly identify 
these Goldstone modes, we will still find bounds on the 
stability. Including more modes will only make the sys- 
tem more unstable. Fig. ^ shows the stability boundary 
for H aa + Q aa and H bb + Q bb in the (A b , A s ) plane. We 
have used 12 texture modes and 6 shape modes in the nu- 
merical calculations. Below are the unstable eigenmodes 
at the stability boundary for the points indicated in Fig. 



(1) aw (0.70,-0.71,0.01,0) 
b« (0.70,-0.71,0.01,0) 

(2) a« (0.04,-0.99,0.01,0) 

b« (0.04,-0.99,0.01,0) (54) 

(3) a« (0,1.0,0) 
bw (0,1.0,0) 

(4) a« (0,1.0,0) 
bw (0,1.0,0) 

We see that when A s is small the first unstable mode 
contains significant contributions from the n = 1 mode 
indicating that the Goldstone modes a = (1,0) and 
b = (1,0) have not been properly identified within the 
numerical approximation. For larger values of A s the first 
unstable modes are pure n = 2 modes and the numerical 
calculation also correctly identifies the n = 1 modes as 
zero-eigenvalue mode. Now looking back at the eigen- 
modes for small A s we note that only the n = 1 and 
11 = 2 modes contribute significantly, suggesting that the 
pure n = 2 mode is always the first mode to go unstable. 



F eff (a, b) = (a\H aa + Q aa \&) + (h\H bb + Q bb \h), (52) 
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6. The stability boundary for (H aa + Qaa) and 
}bb) as found by using 12 texture and 6 shape modes. 



FIG. 
(H bb + 

The two boundaries are apparently the same within the ac- 
curacy of the calculation (« 10%). The unstable eigenmodcs 
associated with the noted points are given in the text, one 
mode each for (H aa + Qaa) and (Hbb + Qbb)- Above these 
points on the boundary the unstable modes appear to be pure 
n = 2 modes. As noted in the text, it seems likely that the 
true unstable modes are always the n — 2 modes. 

Fig. [6] shows us that for any A s > there is a critical 
value of A b below which the system becomes unstable. 
Though not shown in the figure, we remark that as Ab 
becomes even smaller, more and more modes become un- 
stable. This should not be surprising since, if Ab = and 
A s > then all modes are unstable. 

Thus we find that for the free energy (Q) the configu- 
ration consisting of the disk with the preferred exterior 
defect texture [ @ with cp d = and r d = r° d Q] is 
unstable with respect to correlated changes of the shape 
and texture when A s > 0. Furthermore, the absence of 
linear terms in a and b in F arises from a cancellation be- 
tween the bulk and boundary portions of the free energy. 
Hence, changing the structure of either will generically 
generate such linear terms leading to linear shape insta- 
bilities. This helps us understand why either altering the 
form of the boundary energy or accounting for deviations 
from the one-coupling-constant approximation also leads 
to shape instabilities. 



B. The interior defect 

Now we will examine the situation with the interior 
defect. Again, we take D as in (B2h and now we take 



7T 

2' 



(55) 



the preferred texture for a disk as found in section III B| . 
This time, since the defect is inside the disk, it is not 
possible to make an appropriate choice of the P n 's in ( |36| ) 
to change the position of the defect. Nevertheless, we will 
not allow the defect location to vary. This merely has the 
effect of removing the Goldstone modes associated with 
simultaneous translations of the texture and the shape. 

We will again work only to quadratic order in p and 
P and extract an effective free energy F e s(p, K, 7, //, R). 
As with the exterior defect calculation this F c g will con- 
tain no linear terms in p and the quadratic terms will 
not always be positive. We will be able to investigate 
analytically the stability criteria in this case, and we will 
find again that for any A s > and for small enough Ab 
shape instabilities arise. The answer to the question of 
which mode becomes unstable first will be much more 
interesting here. We will find that for the kth harmonic 
mode there is always a value of A s such that this mode 
is the first to go unstable as A b is lowered from 1. It will 
always be assumed that Ab < 1 s ince for larger values 
we already know from section III B that the defect would 



prefer to move away from the origin. This will appear in 
our current analysis as an instability in the n — 1 mode. 

As with the exterior defect case, here again it would 
be necessary to proceed to sixth order in p in the cal- 
culation of -F ff in order to calculate equilibrium shapes. 
Again this task appears formidable and we will only be 
examining the stability of the configurations. 

Proceeding, we have 



F = irK log 



PnP- 



R 2 



neZ-{0} 



neZ-{0} 



PnP- 



R 2 



-2nfiR + tt^lR 

neZ-{0} 
neZ-{0} 
neZ-{0} 

+ 2T7-/R + TT-/R ^2 fa 2-1 ) 
neZ-{0} 



PnP- 

R 2 



(56) 



The absence of any terms linear in P with no p's con- 
firms the extremal nature of the texture. The positive 
definiteness of the quadratic form associated with the 
P's confirms that the texture is in fact a minimum with 
respect to variations that do not change the position of 
the defect. 
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(c) 

FIG. 7. Stability boundaries in the (Ab, A s )-p\ane for a few modes. Note the n = 1 boundary in (a), this mode is the first to 
go unstable as Ab increases beyond 1, instabilities in the other modes for At > 1 are not shown. This instability corresponds 
to the defect wishing to move away from the origin. In all three pictures note the crossing of the stability boundaries. Starting 
in the stable region, as Ab is lowered which mode goes unstable first clearly depends upon the value of A s . For completeness 
we note that the nth stability curve has an asymptote at Ab — 1 — 1/n as one can see in (a) for the n = 2, 3 curves. 



To calculate F e g we simply demand that dF/dPj = 
and solve for Pj(p n ): 



^. = 2n(KijR\j\^ + K\j\P. j R^\ 



yields 



where 



Pa 



\j\ J R^R 



S, = 



K-nR 



K + W R 



(57) 



(58) 



(59) 



Substituting into (pq) we have 



F eS = irK log 



R 



2-KfiR + 2ttjR 



-2?r fiR ^2 



PnP-n 

1 R 2 



0(P 3 



(60) 



where 
Q„ = 



- [Ab - 1 



(A b 



A b + ± 



l + A 



-(n 2 -l) 



(61) 



Note that to 0(p 2 ) F e s is block diagonal in the p„'s and 
diagonal in the corresponding a n 's and b n 's. We can 
examine the stability of the circle simply by looking at 
the sign of Q n . If fi„ is negative then the nth harmonic 
is unstable. 

The stability of the n = 1 mode is independent of A s 
and we can easily see that this mode is unstable when 
Ab > 1. This corresponds precisely to our earlier finding 
( p6| ) that the origin becomes a local maximum for the 
position of the defect when A b > 1. 



For the other modes the boundary between the stable 
and unstable regions in the (Ab, j4 s )-plane is more com- 
plex. For the nth mode this boundary is given by the 
curve Q n = 0, or cquivalently by, 



A° s (n,A b ) 



A,, 



1 - nA h 



1 - A b nA b - n + 1 



(62) 



We are only interested in the region (A s > — 1, < A b < 
1), since physically we require (A b > 0, A s > —1) and we 
have already seen that the n = 1 mode becomes unstable 
when Ab > 1. In Fig. we have plotted A° s (n, A b ) in this 
region for several values of n. Within this piece of the 
(Ab, A s )-plane the nth harmonic mode is stable in the 
region to the right of A° s (n, A b ) and unstable to the left 
of this curve. 

Fixing the physical elastic constants of the system and 
allowing the size of the domain (R) to change is, in the 
(A b , j4 s )-plane, equivalent to fixing A s and allowing A b 
to vary. Thus, a growing domain traces out a horizontal 
line in the (^4^, j4 s )-plane. Starting in the stable region 
(to the right of the stability curves) , as the domain grows 
which stability boundary the horizontal line crosses first 
will clearly depend upon the value of A s . In fact for 



m — 4m + 5m 



4m 



2m 



1 , m 

— <A, < — 



2m - 1 



2m + 1 



(63) 



the stability boundary for the mth harmonic mode will 
be the first such boundary crossed. Thus the domain will 
become unstable with respect to m-fold shape perturba- 
tions first. Of course, as the domain grows further it will 
become unstable with respect to more and more modes. 
Furthermore there is extensive mode coupling at higher 
order in p, and as such the equilibrium shape may not 
exhibit m-fold shape distortions. Nevertheless, this be- 
havior is markedly different from what we found in the 
case of the exterior defect. There it appeared that the 
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n = 2 mode always became unstable first. Interestingly, 
in the extensively investigated dipole model, the n = 2 
mode also appears to be the first mode to become unsta- 
ble. While in both of these cases higher harmonics also 
become unstable as the domain grows larger, and mode 
coupling is present, it is still a possibility that m-fold 
shape distortions may be more easily accounted for by 
the present model with a captured defect. 



V. CONCLUSIONS 

After a brief review of the relevant types of order 
present in domains appearing in the coexistence regions 
of mono-layer films we focused on the effects of the tilt 
order through the simple XY-mode\. We have shown 
the possible existence of a first-order transition in the 
texture for a circular domain from a virtual defect tex- 
ture to a captured defect texture. We then examined the 
shape instabilities for both of these textures. Through a 
slightly different parameterization of the domain we were 
able to reproduce the result of Rudnick and Bruinsma 
|jl6f that the circle is an extremum for the virtual defect 
texture. However, we find that it is not always a mini- 
mum and shape instabilities can arise whenever there is 
the possibility, locally, of a negative effective line tension 
(A a > 0). For the interior defect texture we have found a 
complex stability landscape for the fourier modes associ- 
ated with the domain shape, again arising when A s > 0. 
In both cases larger domains become unstable with re- 
spect to the fourier modes. With the exterior defect tex- 
ture, the n = 2 mode apparently is always the first to 
become unstable. With the interior defect texture which 
mode becomes unstable first depends upon the physical 
parameters of the system. 

While the dipole model can also have shape instabili- 
ties in its fourier modes, it cannot produce chiral domain 
shapes. However, for unequal elastic constants the XY- 
model may yield chiral shapes. Of course, a general sys- 
tem may have important contributions from the dipole 
model and from the XY-mode\. 
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APPENDIX A: SHAPE PARAMETERIZATION 

Here we review some of the properties, including ad- 
vantages and shortcomings, of our choice of parameteri- 
zation for the boundary of the domain D. Recall (p2|) 



r{4>) = Ro 



E 

neZ+ 



(a n cos ncj) + b r , 



(Al) 



Next recall that a domain, D, is said to be star-shaped 
if there exists a point PoeD such that, for all PeD the 
line segment connecting Pq to P is contained in D. The 
boundary of any star-shaped domain, where the origin 
can serve as the point P$eD, can be described by QAl|). 
Correspondingly, a curve given by (Al), if it is also a 
bounded simple closed curve, serves as the boundary of 
a star-shaped domain where the origin can be chosen as 
Po- 

A priori we have no reason to restrict our attention 
to star-shaped domains. Unfortunately our parameter- 
ization does not allow the description of domains that 
are not star-shaped and that have boundaries that are 
simple closed curves (Fig. |8|). Furthermore, the parame- 
terization (Al) admits curves that are not simple closed 
curves and that we would consider to be physically irrel- 
evant (Fig. ||). Thus, if we were to attempt to find highly 
distorted equilibrium domain shapes, we would be well 
advised to use a different parameterization. 






(a) (b) (c) 

FIG. 8. The simple closed curve in (a) is admitted by the 
parameterization ( A_l ) , it bounds a star-shaped domain. In 
(b) this simple closed curve does not bound a star-sh ape d 
domain and is not admitted by the parameterization (Al). 
Finally in (c) we have a closed curve with self-intersections 



which is admitted by (Al) 



In this paper, however, we are only interested in in- 
vestigating the stability of a circular domain. This can 



be accomplished with the parameterization (Al) since 
the circular domain is star-shaped and physically rele- 
vant infinitesimal perturbations will also be star-shaped. 

Another potential difficulty with this parameterization 
is the special role the origin assumes. Namely, we must 
be allowed to choose the origin as our special point Po. 
Thus, for example, a circle of radius 1 centered about 



(2, 0) cannot be described by the parameterization ( Al). 

If we are only concerned with the shape of the domain, 
then, if it is star-shaped, we can always choose the origin 
to be a point in the domain satisfying the properties of 
Po. But now if we are not concerned with the position of 
the domain with respect to the origin of the coordinate 
system, then our parameterization contains redundan- 
cies. Namely, for each point in the domain that satisfies 
the required properties of Po, we can find a parameteriza- 
tion of the form (Al) which produces the correct bound- 
ary. For example, for a disk any interior point could be 
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chosen as the origin. Translating a domain around in gen- 
eral involves changing all of the fourier coefficients and so 
removing this redundancy is not necessarily trivial. How- 
ever, again, being interested only in small perturbations 
of a circular domain allows us to easily deal with this 
redundancy. Starting with a circular boundary centered 
at the origin, 



r(4>) = R, 



(A2) 



one finds that for an infinitesimal transla- 
tion, (e cos ip, e sin-0), the translated boundary is given 
by, 

r{4>) = R + (ecosifj) cos <j> + (e sintp) smcj) + 0(e 2 ). (A3) 

That is, to 0(e) only the n — 1 modes are altered. It is 
for this reason that one often ignores the n = 1 modes 
when using this parameterization. To lowest order these 
modes are merely translations of the domain, and as such 
merely produce annoying redundancies in the parameter- 
ization if one is interested only in the shape of the domain. 
However, in our problems we also have a nontrivial tex- 
ture present which breaks the isotropy of space. Thus the 
position of the domain relative to a fixed origin (provided 
by the texture) is now physically relevant. The redun- 
dancies are now in the simultaneous translations of the 
texture field and the domain. 



APPENDIX B: THE CHOICE OF 77 = 

We would like to comment here on the invariance of 
our results with respect to our choice of boundary energy, 
namely 77 = in (^|). 

Let's consider the two free energies, 

F c = -K I (V$) 2 d 2 x + 77 f m fids 
2 Jd JdD 



-fi f rh • ids + 7 

IdD JdD 



<)D 

ds 



(Bl) 



and 



F a = -K I (V$ fdrx + 77 f rh • hds + 7 f ds. 
2 .Id JdD JdD 



Under the transformation 



$' = *' + tan" 1 H 



(B3) 



we see that F c = F a . Thus, it is trivial to extend our 
results for 77 = to arbitrary values of 77 and fi. 

Interestingly the domain D does not participate in 
the transformation. Thus, although F c contains a chi- 
ral term, indicating that the equilibrium configuration 



(D, rh) should be chiral, we expect D to be achiral. This 
follows from the expectation that the equilibrium config- 
uration for Fa,, (D,rh'), should be achiral since F a con- 
tains no chiral terms. Thus, the chirality only manifests 
itself in the equilibrium configuration through the texture 
rh and not through the shape D. 
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